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Abstract. We describe a novel tool for the quantum characterization of optical 
devices. The experimental setup involves a stable reference state that undergoes 
an unknown quantum transformation and is then revealed by balanced homodyne 
detection. Through tomographic analysis on the homodyne data we are able to 
characterize the signal and to estimate parameters of the interaction, such as the 
loss of an optical component, or the gain of an amplifier. We present experimental 
results for coherent signals, with application to the estimation of losses introduced 
by simple optical components, and show how these results can be extended to the 
characterization of more general optical devices. 



1. Introduction 

Quantum homodyne tomography (QHT) is certainly the most successful technique 
for measuring the quantum state of radiation. It is based on homodyne detection, 
where the signal mode is amplified by the local oscillator. This means that there is 
no need for single-photon resolving photodetectors, whence it is possible to achieve 
quantum efficiency rj approaching the ideal unit value by using linear photodiodes jl[ . 
Moreover, QHT is efficient and statistically reliable, such that it can be used on-line 
with the experiment. Indeed, among other proposed state reconstruction methods, 
QHT is the only one which has been implemented in quantum optical experiments 

SI- 

Possible applications of QHT range from the measurement of photon correlations 
on a sub-picosecond time-scale Q to the characterization of squeezing properties 
j|, |[, photon statistics in parametric fluorescence Q, quantum correlations in down- 
conversion Q and nonclassicality of states In general, the key point is that QHT 
provides information about the quantum state within a chosen sideband, thus allowing 
for a precise spectral characterization of the light beam under investigation. 

In this paper we address QHT as a tool for the quantum characterization of 
optical devices, like the estimation of the coupling constant of an active medium or 
the quantum efficiency of a photodetector. The goal is to link the estimation of 
such parameters with the results from feasible measurement schemes, as homodyne 
detection, and to make the estimation procedure the most efficient. We present our 
experimental results about the reconstruction of the quantum state of coherent signals, 
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together with application to the estimation of the losses introduced by simple optical 
components. Moreover, we show how these preliminary results can be extended to the 
characterization of more general optical devices. 

In the next Section we review some basic elements of quantum tomography, 
whereas in Section ^ we describe the basic requirements needed to implement a 
quantum characterization tool based on tomographic measurements. In Section ^| 
the experimental apparatus is described with some details, and in Section the 
experimental data are analyzed and discussed. Section || closes the paper by discussing 
the possible extensions of the present work. 



2. Quantum Homodyne Tomography 

Quantum tomography of a single-mode radiation field consists of a set of repeated 
measurements of the field-quadrature x^, = \{ae~ % ^ + a^e 1 ^) at different values of the 
reference phase 4>. The expectation value of a generic operator can be obtained by 
averaging a suitable kernel function R[0](x, <j>) as follows 



Tt{qO}= f — f dxp(x. 
Jo n J-oo 



<j>)R[0](x,<f>), (1) 



where p(x, (f) denotes the probability distribution of the outcomes x for the quadrature 
and R[0]{x, 4>) is given by 



1 f°° 

R[0}(x, 4>) = 7 dr Tr {O cos [y/r(x - x*)] } 
4 Jo 



(2) 



Actually, the tomographic kernel R\Q ] (x , 4>) for a given operator O is not unique, 
since a large class of null functions |9|, [10| F(x, (j>) exists that have zero tomographic 
average for arbitrary state. This degree of freedom can be exploited to adapt the 
kernel to data and achieve an optimized determination of the quantity of interest. For 
example, quantities like the photon number, the field amplitude and any normally 
ordered moment can be the tomographically estimated by averaging the following 
kernel 



M-l M-l 



K[0]{x,4>)=R[0]{x 1 4>)+ > u k F k (x,<t>)+ > u* k F*{x,4>) , (3) 



k=0 k=0 



where the kernels R[0] for the moments are given by [[Tl| 

R[a^ n a m ] (x; 0) = e ^(™-«)0 J^+^M ( 4 ) 

H n (x) being the Hermite polynomial of order n, and the null functions F k are 
expressed as Fk(x,<f>) = x k e^ k+2 ^ , k — 0,1,... The coefficients fj,k are obtained 
by minimizing the rms error for the given kernel on the given homodyne sample. A 
similar approach can be applied to optimize the reconstruction of the matrix elements 
Qmn = ( m |f?|^)j thus achieving an effective quantum state characterization. 



3. QHT as a tool to estimate parameters 

The state reconstruction method provided by QHT is effective and reliable, such that 
it can be exploited to build a tool to characterize quantum devices. The general scheme 
for such a tool should be as follows. First, we need a stable source of quantum states, 
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i.e. a source able to provide repeated preparations of a reference signal. The signal 
can be characterized by QHT, and then employed as input of a given device, which we 
want to characterize by the estimation of some relevant parameters. At the output, 
the transformed state can analyzed by QHT, such to characterize the input-output 
relations of the device. 

In order to implement this kind of scheme two basic requirements should be 
satisfied: i) we need a stable source for the reference signal, and ii) an effective data 
processing for the tomographic samples should be devised, in order to minimize the 
number of measureements. 

The first point can be satisfied by considering Gaussian signals, like coherent 
or squeezed states. Indeed, quantum signals that are most likely to be reliably 
generated in a lab are Gaussian states. The most general Gaussian state can be 
written as g = D(n) S{r)v (r) (fi) , where v denotes a thermal state v = 
(n t h + i)^ 1 [nth/(nth + l)] a a , S(r) = cxp[r(a 2 — a^ 2 )/2] the squeezing operator and 
Z?(/i) = exp(/ia^ — /i*a) the displacement operator. However, thermal excitations 
can be neglected at optical frequencies, such that we may generally consider v as the 
vacuum state. The homodyne distribution of the state g at phase (f> with respect to 
the local oscillator is Gaussian and, remarkably, such Gaussian character is not altered 
by many transformations induced by optical devices, such as the loss of a component, 
the gain of an amplifier or the quantum efficiency of a detector. In this paper we 
consider the reference signal excited in a coherent state. More general signals will be 
considered elsewhere. 

The need of an effective data processing lead to consider either adaptive or 
maximum-likelihood (ML) procedures on the tomographic data. In Section ^ we apply 
adaptive tomography for the estimation of losses induced by optical filters. Here, we 
illustrate the use of ML procedure to the characterization of a general (active or 
passive) optical media, which we plan to perform experimentally in the near future. 

Let us start by reviewing the ML approach. Let p(x\\) the probability density 
of a random variable x, conditioned to the value of the parameter A. The analytical 
form of p is known, but the true value of the parameter A is unknown, and should 
be estimated from the result of a measurement of x. Let X\,X2, ...,Xjsr be a random 
sample of size N. The joint probability density of the independent random variable 
Xi,X2, :.,Xi\f (the global probability of the sample) is given by 

£(A)=nf =1 p(x fc |A), (5) 

and is called the likelihood function of the given data sample. The ML estimator of 
the parameter A is defined as the quantity Ami that maximizes £(A) for variations of 
A. Since the likelihood is positive this is equivalent to maximize 

N 

L(X)^\ogC(X)=J2^gp(x k \X) (6) 

k=l 

which is the so-called log-likelihood function. 

Let us consider a generic optical media: the propagation of a signal is governed, 
for negligible saturation effects, by the master equation 

g = GiL[a]g + G 2 L[a ] ]g , (7) 

where g is the density matrix describing the quantum state of the signal mode a and 
L[0] denotes the Lindblad superoperator L[0]A = O AO^ - \O^OA- \AO ] 0. If 
we model the propagation as the interaction of a traveling wave single-mode a with a 
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system of N identical two-level atoms, then the absorption G\ = "fNi and amplification 
G2 = 7^2 parameters are related to the number Ni and N 2 of atoms in the lower and 
upper level respectively. The quantity 7 is a rate of the order of the atomic linewidth 
p2[ , and the propagation gain (or deamplification) is given by Q = exp[(G*2 — G\) t\. 
A medium described by the master equation (^) represents a kind of phase-insensitive 
optical device, such that the parameters G\ and G2 can be estimated starting from 
random phase tomographic data. According to (|7|), the homodyne distribution of a 
coherent signal with initial amplitude ao is given, after the propagation, by 

= ^ 1 +gV2) ^{-^V2 [*-«^)] 2 } W 

with g = e~ Qt , 2Q = {G x - G 2 ), and S 2 = (d + G 2 )(l - g 2 )/AQ (for non-unit 
quantum efficiency rj < 1, 6 2 + g 2 /2 — ► 5 2 + g 2 /2 + (1 — 77) / 277) . By ML estimation 
on homodyne data we may reconstruct the parameters G\ and G2. The resulting 
method has proven efficient by numerical simulations JL3| , and provides a precise 
determination of the absorption and amplification parameters of the master equation 
using small homodyne data sample. Notice that no advantage in using squeezed states 
should be expected, because of the phase-insensitive character of the device. 

4. Experiment 

A schematic of the experiment is presented in Fig. [l| The principal radiation source 
is provided by a monolithic Nd:YAG laser (« 50 mW @1064 nm; Lightwave, model 
142). The laser has a linewidth of less than 10 kHz/ms with a frequency jitter of less 
than 300 kHz/s, while its intensity spectrum is shot-noise limited above 2.5 MHz. 
The laser emits a linearly polarized beam in a TEM00 mode. 



FR HWP, HWPj 




Figure 1. Schematic of the experimental set-up. A Nd:YAG laser beam is 
divided into two beam, one acts as the homodyne local oscillator, the other 
represent the signal beam. The signal is modulated at frequency Q. with a defined 
modulation depth to control the average photon number in the generated coherent 
state. The tomographic data are provided by an homodyne detector whose 
differece photocurrent is demodulated and then acquired by a digital oscilloscope 
(Tektronix TDS 520D). 



The laser is protected from back reflection by a Faraday rotator (FR) whose 
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output polarizing beam splitter can be adjusted to obtain an output beam of variable 
intensity. The beam outing the isolator, of w 2.5 mW, is then split into two part 
of variable relative intensity by a combination of a halfwave plate (HWPi) and a 
polarizing beam splitter (PBS, p-transmitted, s-reflected) . The strongest part, is 
directly sent toward the homodyne beam-splitter (BS) where it acts as the local 
oscillator beam. One of the mirror in the local oscillator path is piezo mounted 
to obtain a variable phase difference between the two beam. The remaining part, 
typically less than 200 /«W, is the homodyne signal. The optical paths traveled by 
the local oscillator and the signal beams arc carefully adjusted to obtain a visibility 
typically above 75% measured at one of the homodyne output port. The signal beam 
is modulated, by means of a phase electro-optic modulator (EOM, Linos Photonics 
PM0202), at 4MHz, and a halfwave plate (HWP 2 , HWP 3 ) is mounted in each path 
to carefully match the polarization state at the homodyne input. 

The basic property of the homodyne detector (described into details in Ref. p5[ ) 
is a narrow-band detection of the field fluctuations around 4MHz. The detector is 
composed by a 50-^50 beam splitter (BS), two amplified photodiodes (PD1, PD2), 
and a power combiner. The difference photocurrent is demodulated at 4MHz by 
means of an electrical mixer. In this way the detection occurs outside any technical 
noise and, more important, in a spectral region where the laser does not carry excess 
noise. 

The phase modulation added to the signal beam move a certain number of 
photons, proportional to the square of the modulation depth, from the carrier optical 
frequency lu to the side bands at uj ± SI so generating two few photons coherent state, 
with engineered average photon number, at frequencies u> ± Q. The sum sideband 
modes is then detected as a controlled perturbation attached to the signal beam J3[. 
The demodulated current is acquired by a digital oscilloscope (Tektronix TDS 520D) 
with 8 bit resolution and record length of 250000 points per run. The acquisition is 
triggered by a triangular shaped waveform applied to the PZT mounted on the local 
oscillator path. The piezo ramp is adjusted to obtain a 2n phase variation between 
the local oscillator and the signal beam in an acquisition window. 

The homodyne data, to be used for tomographic reconstruction of the radiation 
state, have been calibrated according to the noise of the vacuum state. This is obtained 
by acquiring a set of data leaving the signal beam undisturbed while scanning the local 
oscillator phase. It is important to note that in case of the vacuum state no role is 
played by the visibility at the homodyne beam-splitter. 

5. Data Analysis 

Our tomographic samples consist of N homodyne data {xj, 4>j}j=i,...,N with phases 4>j 
equally spaced with respect to the local oscillator. Since the piezo ramp is active during 
the whole acquisition time, we have a single value xj for any phase cf>j. From calibrated 
data we first reconstruct the quantum state of the homodyne signal. According to the 
experimental setup described in the previous section we expect a coherent signal with 
nominal amplitude that can be adjusted by varying the modulation depth of the 
optical mixer. However, since we do not compensate for the quantum efficiency of 
photodiodes in the homodyne detector {-q ~ 90%) we expect to reveal coherent signals 
with reduced amplitude with respect to actual one. In addition, the amplitude is 
furtherly reduced by the non-maximum visibility (ranging from 75% to 85%) at the 
homodyne beam-splitter. 
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In Fig. ^| we show a typical reconstruction, together with the reconstruction of 
the vacuum state used for calibration. For both states, we report the raw data, the 
photon number distribution, i.e. the diagonal elements g nn = (n|g|n) of the density 
matrix in the Fock representation, and a contour plot of the Wigner function. The 
matrix elements are obtained by sampling the corresponding kernel functions 



R[ \n)(n + k\ ] (x, <f>)=2 exp {-ik<f>} \J 2 k n\(n + k)\ f nk (x) , 



where 
fnk(x) 



ELo + l + *A V2; -2, 2 ) * even 



e;u { - y ^;^;: k T )m ^-^H-i - * A ^ **) * odd 



l\{n-l)\{l+k)\ 

and $(a, b;x) denotes a confluent hypergeometric function. The tomographic 
determination of the matrix elements is given by the averages 

Qnk = R[ \n)(k\ ] = jjY^R[ \n)(k\ Kay, <fo) , (9) 

3 

whereas the corresponding confidence intervals are given (for diagonal elements) by 
fiQnn — Ag/y/N, Ag being the rms deviation of the kernel R over data (for off- 
diagonal elements the confidence intervals are evaluated for the real and imaginary 
part separately). 

In order to see the quantum state as a whole, we also report the reconstruction 
of the Wigner function of the field, which is defined as follows 

W{z)^^[ Q D(2z)(-f a ) , (10) 

and can be expressed in terms of the matrix elements as 

oo oo 

W(z) = ReJ2 ^ E A ( n ' d ' l*l a K«+d (11) 



d=Q n=0 



where 



A(n,d; M 2 ) = (-)"2(2 - ^ )|2z|^^- i -^ ye -^|-Lf l (|2z| 2 ) , (12) 

and L^{x) denotes the Laguerre polynomials. Of course, the series in Eq. ( fTl| ) should 
be truncated at some point, and therefore the Wigner function can be reconstructed 
only at some finite resolution. 

Once the coherence of the signal has been established we may use QHT to 
estimate the loss imposed by a passive optical component like an optical filter. The 
procedure may be outlined as follows. We first estimate the initial mean photon 
number ho = |c*o| 2 of the signal beam, and then the same quantity inserting an 
optical neutral density filter in the signal path. If T is the loss parameter, then the 
coherent amplitude is reduced to ar — aoe _r , and the intensity to t%y — hae" 2r . 

The estimation of the mean photon number can be performed adaptively on data 
by taking the average of the kernel 

K[<Ja]{x) = 2x 2 -^+ ^ + fi*e- 12 * , (13) 

where fj, is a parameter to be determined in order to minimize fluctuations. As proved 
in Ref. [pi \i = — 1/2 (a t2 ), which itself can be obtained from homodyne data. In 
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Figure 2. Reconstruction of the quantum state of the signal, and of the vacuum 
state used for calibration. For both states, from left to right, we report the raw 
data, a histogram of the photon number distribution, and a contour plot of the 
Wigner function. Reconstruction has been performed by a sample of TV = 242250 
homodyne data. The coherent signal has an estimated average photon number 
equal to (at a) = 8.4. The solid line denotes the theoretical photon distribution 
of a coherent state with that number of photons. Statistical errors on matrix 
elements are about 2% . The slight phase asymmetry in the Wigner distribution 
corresponds to a value of about the 2% of the maximum. 

practice, one uses the data sample twice: first to evaluate /i, then to obtain the 
estimate for the mean photon number. 

In Fig. [j] the tomographic determinations of rip are compared with the expected 
values for three set of experiments, corresponding to three different initial amplitudes. 
The expected values are given by np = noe~ 2r V, where T is the value obtained 
by comparing the signal dc currents Iq and Jr at the homodyne photodiodes and 
V = Vr /Vo is the relative visibility. The solid line in Fig. || denotes these values. The 
line is not continuos due to variations of visibility. As it is apparent from the plot the 
estimation is reliable in the whole range of values we could explore. It is worth noticing 
that the present estimation is absolute, i.e. it does not depends on the knowledge of 
the initial amplitude, and it is robust, since it may performed independently on the 
quantum efficiency of the photodiodes employed for the homodyne detector. 

6. Conclusions and outlook 

In this paper we addressed QHT as a tool for the characterization of quantum optical 
devices. We carried out the quantum state reconstruction of coherent signals, and 
show how QHT can be used to reliably estimate the loss imposed by an optical filter. 
We also show that the estimation procedure can be extended to the characterization 
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Figure 3. Estimation of the mean photon 
number of a coherent signal as a function of 
the loss imposed by an optical filter. Three set 
of experiments, corresponding to three different 
initial amplitudes are reported. Open circles are 
the tomographic determinations, whereas solid 
line denotes the expected values, as follow from 
nominal values of loss and visibility at the 
homodyne. Statistical errors are within the circles. 
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of general (active or passive) optical devices, which we plan to perform experimentally 
in the near future. We also plan to extend our analysis to squeezed signals since it 
has been proved fl4|| that squeezing improves precision in the homodyne estimation of 
relevant parameters such the phase-shift or the quantum efficiency of a photodetector. 
In particular, our aim is to reproduce the described measurement scheme at the output 
of an OPO cavity jl6| driven below threshold to generate vacuum squeezed radiation. 
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